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Do diffusive non-Markovian stochastic Schrodinger equations (SSEs) for open quantum systems 
have a physical interpretation? In a recent paper [Phys. Rev. A 66, 012108 (2002)] we investigated 
this question using the orthodox interpretation of quantum mechanics. We found that the solution 
of a non-Markovian SSE represents the state the system would be in at that time if a measurement 
was performed on the environment at that time, and yielded a particular result. However, the 
linking of solutions at different times to make a trajectory is, we concluded, a fiction. In this paper 
we investigate this question using the modal (hidden variable) interpretation of quantum mechanics. 
We find that the noise function z(t) appearing in the non-Markovian SSE can be interpreted as a 
hidden variable for the environment. That is, some chosen property (beable) of the environment has 
a definite value z(t) even in the absence of measurement on the environment. The non-Markovian 
SSE gives the evolution of the state of the system "conditioned" on this environment hidden variable. 
We present the theory for diffusive non-Markovian SSEs that have as their Markovian limit SSEs 
corresponding to homodyne and heterodyne detection, as well as one which has no Markovian limit. 
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I. INTRODUCTION 

In nature it is very unlikely to find a system existing 
in isolation; usually it is immersed in an environment 
(or bath). In quantum mechanics we label this type of 
system an open quantum system [1]. To determine the 
evolution we must solve the Schrodinger equation 



-H uni (t)\*(t)}, 



(1.1) 



where and H un [(t) are the quantum state and 

Hamiltonian for the complete universe (system and 
bath). That is, [^(i)) belongs to the Hilbert space 



'H-a 



Wbath- Due to the large Hilbert space of 



the bath (7ibath) it is convenient to describe the system 
by its reduced state. This is defined as 



fred(t)=nWh[|*(f)>(*(f)|], 



(1.2) 



and operates only in the Ti. sys . 

It has been shown [2, 3] by a projection-operator 
method that we can write a general master equation for 
the reduced state as 



dtPred{t) = --[H(t),p Tcd (t)] + 
n 



K(t,s)pred(s)ds, 



(1.3) 

where H(t) operates only in H. sys and K,{t, s) is the "mem- 
ory time" superoperator. It operates on p le d(s) and rep- 
resents how the bath has affected the system in the past. 
The problem with this equation is that in general K,(t, s) 
cannot be explicitly evaluated. 



In recent papers [4, 5, 6, 7, 8, 9, 10, 11, 12, 13] non- 
Markovian stochastic Schrodinger equations (SSEs) have 
been proposed which allow an alternative procedure for 
solving the reduced state. A non-Markovian SSE is a 
stochastic equation for the system state \ip z (t)), condi- 
tioned on some noise function z{t). We consider only 
continuous SSEs, although discontinuous ones have also 
been proposed [14] . The SSE has the property that when 
the projector for \tp z (t)} is averaged over all the possible 
z(i) one obtains /o re d(t). That is, 



Pied (t) = vu z (t))(Mt)\l 



(1.4) 
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where E[...] denotes an ensemble average over all possible 
z(t)'s. 

When using non-Markovian SSEs to solve the reduced 
state it turns out that in general we cannot explicitly 
evaluate \ip z (t)). However, as shown in Refs. [15, 16] 
we can use perturbative techniques to find approximate 
solutions. Here we are interested not in how to solve the 
non-Markovian SSEs, but in how to interpret them. 

When Markovian SSEs (the Markovian limit of non- 
Markovian SSEs) were introduced in quantum optics 
their first interpretation was as a numerical tool [17] . An- 
other interpretation was that Markovian SSEs represent 
objective (that is, independent of any observer) trajec- 
tories for the system [18, 19, 20]. In this interpretation 
the SSE is seen as a (stochastic and nonlinear) modifi- 
cation of the Schrodinger equation, modelling how state 
reduction (collapse of the wave function) occurs in open 
quantum systems. However, in recent times it has been 
generally accepted that Markovian SSEs are evolution 
equations for the system state conditioned on continuous 
measurement of the bath [1, 21, 22, 23]. 

A Markovian bath is one for which the bath corre- 
lation time is much less than the decoherence time of 
the system. For such a system one can envisage mak- 
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ing repeated measurements on the bath on a time scale 
infinitesimal compared to the system decoherence time 
but large compared to the bath correlation time. Thus 
these measurements do not lead to any back-action of 
the system from the bath. But the measurement of the 
bath does yield information about the system. This can 
be described as the action of a measurement operator 
(as appears in generalized quantum measurement theory 
[24] ) for an infinitesimal time interval [23] , and keeps the 
system in a pure state. It can be reexpressed as a SSE in 
which the stochastic variable z(t) is related to the results 
of the continuous measurements. The stochastic trajec- 
tory of this conditioned system state generated by the 
Markovian SSE is often referred to as a quantum trajec- 
tory [1]. Different detection procedures (measurements 
on the bath) result in different unravelings. A few com- 
mon examples arc direct detection [1], homodyne [1, 21], 
and heterodyne detection [22]. 

In the light of the quantum trajectory interpretation 
of Markovian SSEs, we return to the interpretation of 
non-Markovian SSEs. Obviously one interpretation is 
that they are simply a numerical tool used to generate 
the reduced state. However, after the success in find- 
ing a physical interpretation for Markovian SSEs, it is 
natural to seek something beyond this trivial interpre- 
tation. Moreover, we have previously shown [10] that 
there are different non-Markovian unravelings, and that 
these correspond to different measurement schemes (ho- 
modyne and heterodyne) in the Markovian limit. Thus 
it is natural to seek an interpretation of non-Markovian 
SSEs beyond that of being just numerical tools [25]. 

In Rcf. [10] we came to the conclusion that under the 
orthodox interpretation [24] of quantum mechanics, the 
solution of a non-Markovian SSE at time t is the state 
the system will be in, if at that time a measurement was 
performed on the bath and yielded a result z. Thus the 
non-Markovian SSE under this view has no interpreta- 
tion; it is just a stencil used to calculate the conditioned 
system state at a particular time t. In other words, the 
linking of \ip z (t)) (or z(t)) with itself at times less then t 
turns out to be a convenient fiction. 

Unlike the Markovian case, it is not possible to de- 
rive the SSE by continuously measuring the bath because 
a non-Markovian bath has a non-negligible correlation 
time. Thus if a measurement is made at time t, collaps- 
ing the bath state at that time, this will affect the state of 
the bath interacting with the system in the future. That 
is, the measurement on the bath will cause a backaction 
on the system and hence the average evolution for a sys- 
tem state conditioned on continuous measurement of the 
bath will not be that of Eq. (1.2). In fact, a continuous 
measurement on a bath with a nonzero correlation time 
will lead to a quantum Zeno effect, radically altering the 
average evolution of the system. If, on the other hand, 
the measurement on the bath is not done continuously, 
then the system state will not remain in a pure state as 
the system and bath will become entangled in the time 
between measurements on the bath. Thus it seems safe to 



conclude that there is no continuous measurement inter- 
pretation, and in fact we are forced to accept that SSEs 
are only a numerical tool which could be used to deter- 
mine the system state at a particular time conditioned 
on a particular measurement result z. 

Since orthodox quantum mechanics fails to give a sat- 
isfying interpretation for non-Markovian SSEs, in this 
paper we turn to a nonorthodox approach: the modal in- 
terpretation of quantum mechanics [26, 27, 28, 29, 30, 31, 
32]. This interpretation, unlike the orthodox interpreta- 
tion, has as its basic goal to keep reality intact. That 
is, the values of some observables (the hidden variables) 
really exist before we measure them. Because the observ- 
ables have an objective reality from now on we will refer 
to them as properties or beables (after Bell [26]). Just as 
in the orthodox theory, where it is impossible to simulta- 
neously measure all observables, in the modal theory it 
is impossible to give all observables property status The 
best- known example of such an interpretation is Bohmian 
mechanics for particles [33] in which position is the pre- 
ferred observable (property). 

We expect a modal interpretation to be applicable to 
non-Markovian SSEs because we can use it to assign defi- 
nite properties to the bath, as occurs in the orthodox the- 
ory when the bath is measured, without invoking such a 
measurement. In this way we avoid the backaction prob- 
lem which arose in the orthodox theory. While the bath 
is ascribed definite properties, the system is described as 
a purely quantum system. But, because of the entangle- 
ment between the system and the bath, we can define a 
system state associated with (or "conditioned" on) a par- 
ticular value for the bath property. If the bath properties 
are described by rank-one projective measures on 7ibath, 
then the conditioned system state will be pure. Averag- 
ing over the conditioned system state would reproduce 
the non-Markovian reduced state matrix p Te d(t), just as 
in Eq. (1.4). 

Since the value of the bath hidden variables change 
in time, the conditioned system state will evolve in time 
also, and this time there is a meaningful relation between 
the conditioned system state at different times. A par- 
ticular set of bath properties amounts to a particular 
decomposition (of the unit operator on 7ibath), and leads 
to a particular unraveling of the non-Markovian master 
equation (1.3). We will use these terms interchangeably. 

In this paper we show that for a suitable choice of bath 
properties we can reproduce all of the non-Markovian 
SSEs discussed above, as well as one which has not been 
previously considered. The noise function z(t) appear- 
ing in the non-Markovian SSE is simply a function of the 
values of the bath hidden variables. The system state 
\i[> z (t)) is the system state conditioned on the bath prop- 
erties having the values giving z(t). Moreover, the sys- 
tem state guides the values of the hidden variables of the 
bath. 

The format of this paper is as follows: In Sec. II a 
summary of the modal interpretation of quantum me- 
chanics is presented. In Sec. Ill we give the microscopic 
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model for non-Markovian dynamics for an open quan- 
tum system that underlies all of our work. In Sec. IV we 
show how modal dynamics can be applied to such mod- 
els, and derive the hidden variable interpretation of non- 
Markovian SSEs for three different unravclings. These 
correspond to objective values for the position of each 
bath harmonic oscillators (as in the Bohmian interpreta- 
tion), for the quadrature of the collective bath field, and 
for the coherent amplitude of the collective bath field. 
The first of these has no Markovian limit, while the lat- 
ter two have homodyne and heterodyne detection as their 
Markovian limit. Lastly in Sec. V we conclude with a 
discussion and directions for future work. 



II. MODAL INTERPRETATION OF QUANTUM 
MECHANICS 

In this section we give a brief overview of the modal in- 
terpretation of quantum mechanics; for a more detailed 
description see Refs. [26, 27, 28, 29, 30, 31, 32]. The 
basic idea of this view of quantum mechanics is that cer- 
tain observables have a definite reality independent of 
measurement, whereas in the orthodox interpretation the 
measurement defines the reality of the observable. To ex- 
plain this we consider measurement of the observable Z . 
The notation we use to represent an observable is 

Z = {{z ni % n )}. (2.1) 

That is, Z is represented by a set of paired elements. 
Each pair represents the value z n and the corresponding 
projector % n . The projectors are orthogonal and form a 
decomposition of unity: 

5>„ = 1 ( 2 - 2 ) 

n 

In Eq. (2.1) for simplicity we have only considered the 
case where Z is time independent. We can in this case 
also define an operator Z which is equivalent to Eq. (2.1), 
by 

Z = ^Z„7T n . (2.3) 

n 

In the orthodox interpretation Z has a definite value 
z„ if and only if \^(t)) = \^„(t)) (an eigenstate of Z). 
However, in general |\&(t)) = Cn(t)\^f n (t)), which im- 
plies that Z has all values contained in the sum; it is not 
well defined. Upon measurement, by the introduction of 
a reduction equation [|*(t)) — > |^Ti(i))] the value of Z 
becomes z n with probability |c n | 2 . 

In the modal interpretation we choose one projective 
measure, {7r ra }, as the preferred measure. This then de- 
termines which observables can be given property status. 
To explain why the property takes the value z n we intro- 
duce an extra quantum state, the property state, 



where N is a normalization constant. This state propa- 
gates in time along with \^(t)), except it is stochastic in 
nature (jumps between different n). It is interpreted as 
the actual state of the universe, by the eigenstate eigen- 
value it has a definite value (z n ) for Z. The stochastic 
dynamics (rates at which it jumps between different n) is 
determined by | and so in this interpretation 
is called the guiding state. 

For the purposes of this paper Eq. (2.3) is actually not 
general enough. In Ref. [32] we showed that this theory 
can be extended to positive operator measures (POMs) 
[34], that is 

Z = {(z n ,F n )}, (2.5) 

where F n is an effect (or POM element) with ^ n F n = 1. 
In Eq. (2.5), z n is the value of the effect, which could be 
a real number, a complex number, a string of numbers, 
or even a statement (yes/no). 

In Ref. [32] we showed that by implementing 
Naimark's theorem, modal dynamics can be extended 
to include POMs. Naimark's theorem says that if wc 
enlarge the Hilbcrt space of the universe from H. unl to 

= H U ni®W a ux) we can define a projector n n in K such 
that 

(*(t)|# n |*(t)) = (*(t)|^|n n |0)|*(t)> > (2.6) 

for all \^{t)) G H U ni and for all possible n. \<fi}(4>\ is called 
the Naimark projection of JC onto 7i un i- To work out the 
set {n„(t)} it is necessary to introduce another projector 
n A r + i(i), such that 

N+l 

ft n (t) = iuni+aux, (2-7) 

a 

and 

n„(t)n m (t) = u n (t)s 

is satisfied for n, m,=l, N+l. The set of projectors 
in this enlarged Hilbert space is called the Naimark ex- 
tension of F n (t) [34] . Worked examples of this are shown 
in [34]. In this enlarged space the observable defined by 
Eq. (2.5) becomes the property 

Z = {(z„,n„)}, (2.9) 

or equivalently 

N+l 

Z=^Z n tl n . (2.10) 

Here {n„} is the preferred projective measure in IC. The 
guiding state becomes 



\* n (t))=w n \*(t))/y/N, 



(2.4) 



i*(t)> = m*)) ® (2.H) 
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where is still the solution to the Schrodingcr equa- 

tion (1.1). The property state becomes 

|$ n (t)) = fl n mt))/VN. (2.12) 

That is, when describing POMs the universe is described 
by a property state in the enlarged Hilbert space, which 
does not factorize into a universe state and an auxiliary 
state. This can be regarded merely as a mathematical 
construction to give reality to POMs. Note that the value 
of zjv+i is irrelevant as IIjv+i projects into the null space 
of |<&(£)). That is this value will occur with zero proba- 
bility. 

The modal dynamics [the stochastic evolution of the 
states in Eq. (2.4) or Eq. (2.12)] is found using the 
method originally proposed by Bell [26] and generalized 
in Refs. [29, 30, 32] to include time dependent projectors 
and POMs. Define P n (t) as the probability that the sys- 
tem is in the n th state at time t. Assuming a Markovian 
process, by which we mean that the probability of being 
in state m at time t + dt only depends on the state at 
time t, we can write a master equation for P n (t) as 

d t Pn(t) = Y^l T nm(t)P m (t) - T mn P n (t)], (2.13) 
m 

where T nm are transition rates. For n = m, T nn (which is 
negative) is a measure of the rate at which state n losses 
probability. 

Defining a probability current J nm (t) as 



(2.14) 



results in J n m(t) = —Jmn(t) and allows us to rewrite the 
probability master equation as 



d t Pn(t) 



(2.15) 



Given J nm {t) and P n (t), there are many possible transi- 
tion rates satisfying Eq. (2.15). One solution, chosen by 
Bell [26] is as follows. 
For J nm {t) < 0, 

T nm {t) = 0, (2.16) 
T mn (t) = -J nm (t)/P n (t), (2.17) 

and for J nm (t) > 

T nm (t) = J nm (t)/P m (t), (2.18) 
T mn (t) = 0. (2.19) 

This is only one of the infinitely many solutions. These 
are found by adding an extra term T° m (t) to T nm (t), 
where T® m (t) is constrained only by 



T° m (t)P m (t)-T^ n (t)P n (t) = 0. 



(2.20) 



In Rcf. [32] we showed that one possible solution for 

Jnmij^) is 



J nm (t) = 2Im{($(i)|n„[ J ff uni (i)® i £ 

xfi m mt))}/h. 



Note that this is only one of infinitely many possible cur- 
rents, as we can add any current J n ° m (t) to J nm (t) which 
satisfies J n ° m = 0, to give a valid probability current. 
For the purposes of this paper we only consider the sim- 
ple solutions [not containing the extra T^ m (t) and Jn m (t) 
terms] . 

The above dynamics only describe a discrete decompo- 
sition. For non-Markovian SSEs we must consider con- 
tinuous decompositions. In Ref. [32] we showed that 
provided the Hamiltonian is at most quadratic in the 
conjugate variable to the chosen hidden variable, then 
the above dynamics reduces to a deterministic theory. 
The only stochasticity is due to the random initial con- 
ditions of the hidden variable. That is, it is similar to 
Bohmian mechanics [33]. In fact if the preferred pro- 
jective measure is chosen to be the position projective 
measure (jt x = |a;)(a;|) then Bohmian mechanics is the 
continuum limit. 

We define continuous decompositions by the preferred 
projective measure {II{ gfc } = !{%}) ({<Zfc}|}, which in turn 
defines {cj k } as the operators for the hidden variables. 
Here we have introduced a notation for the set of hidden 
variables ({<Zfe}) whose relevance will become important 
in the following sections. The decomposition of unity is 



\Uk})(Uk}\~[[d qk = i, 



(2.22) 



For this set of hidden variables we will see later that it is 
also useful to define a single property Z(t) as 



Z(t.) = {(z({ qk },t),\{ qk })({ qk }\)}, 



(2.23) 



where z({ qk },t) is a function of the values of the hidden 
variables { qk } and t. 

As in Bohmian mechanics the evolution of each value 
of the hidden hidden variable qk (t) (or the corresponding 
property state) can be found by the differential equation 



d t qk{t) =v k ({qj},t)\ qk = qk ( t ), 



(2.24) 



where v k ({<jj}, t) is the velocity field. This is defined as 
[32] 

Vk{{q ^ t} - <*(t)i{*}><te}i*(t)> ■ (2 - 25) 

Here v k (t) is the k th component of the velocity operator 
and is defined as 



v k (t) = -^[ qk ,H nni (t) O l a 



(2.26) 



III. 



NON-MARKOVIAN DYNAMICS FOR 
OPEN QUANTUM SYSTEMS 



(2.21) 



The aim of this section is to outline the underlying 
dynamics used to generate non-Markovian SSEs. Firstly 
we assume that there are only two systems, a system 
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of interest and a bath. That is, the Hamiltonian of the 
universe is 

H Bys <g> i + i <g> i?bath + V. (3.1) 

The system Hamiltonian H sys is split into two terms, 
these being Hq and H. The bath is modelled by a col- 
lection of one-dimensional harmonic oscillators. In terms 
of the bath annihilation and creation operators, hk and 
a k , the Hamiltonian for the bath is 

-ffbath = fiuj k a\a k . (3.2) 

k 

The interaction Hamiltonian V we assume is linear in 
the bath amplitude, and has the form 

V = ihY,[L9*A-ttg k a k ], (3.3) 

fe 

where L is the system lowering operator and gk is the 
coupling strength of the fc th mode to the system. With- 
out loss of generality we can take the g k & to be real, 
absorbing any phases in the definitions of the bath oper- 
ators. 

For calculation purposes we define an interaction frame 
such that the fast dynamics placed on the state by the 
Hamiltonians Hq and -ffbath is moved to the operators. 
The unitary evolution operator for this transformations 
is 

U o (t,0) = e-^n^i+i^bathXt-o)/^ /g ^\ 

Thus the combined state in the interaction frame is define 
as 

|*(t)> 0)|*(t)sch>, (3.5) 

and an arbitrary operator A becomes 

A ixA (t) = U^t,0)Au o {t,0). (3.6) 
This allows us to write the Schrodinger equation as 

d t \^(t)) = -^H uni (t)\^(t)), (3.7) 

where H uni (t) = H int (t) + Vmt(t)- Here H int (t) refers to 
H{t) in the interaction picture and the form of the latter 
is 

MntW = ih^lLgke^al - Dgke-^ak], (3.8) 
fe 

where VLk = cu k — SI. Here we have finally restricted the 
form of Hq to be such that L in the interaction picture 
simply rotates in the complex plane at frequency SI. That 
is L int (t) = Le- mt . 



IV. NON-MARKOVIAN STOCHASTIC 
SCHRODINGER EQUATIONS DERIVED 

In this section we show that diffusive non-Markovian 
SSEs have an interpretation under the modal interpre- 
tation of quantum mechanics. To do this we choose a 
decomposition such that the preferred projectors have 
the form 

*{q k ] = |{<7fc})({<7fc}|bath® Isys- (4.1) 

This means the bath is given definite properties, while 
the system is treated as a purely quantum system, which 
nevertheless influences the bath values via the coupling 
Hamiltonian. The different unravelings correspond to dif- 
ferent choices of \{q k }) ({<Zfc}|bath- For the overcomplctc 
unravelings, like the coherent state unraveling [5, 10] we 
have to use a POM [32]. This means we have to use the 
enlarged Hilbert space preferred projector Tl{ qk j, which 
in general form is 

Hfe} = |{2fc})({<Zfc}|bath+aux ® Isys- (4.2) 

With these projectors, the property states \^> n (t)) 
[similarly for |$„(£))] can be factorized as 

!*»(*)> = \Uk})\tp{ qk }{t)), (4.3) 

where \4>{q k } (t)) is called the conditioned system state. 
It receives this name because it lives entirely in Ti. sys 
and is conditioned on the bath values {qk}- The form 
of IVW(*)) is 

|VW*)> = <{<fc}|*(*)>/VJV, (4.4) 

where the normalization constant is defined as 

N={*(t)\{q k })({ qk }\*(t)). (4.5) 

This is the state of the system conditioned on the bath 
hidden variables having values {q k }. 

For an actual trajectory (in the sense of Rcf. [10]), the 
bath values {q k (t)} are time dependent. This state be- 
comes \tp{q k (t)} (t)) and represents the state of the system 
conditioned on the bath having this trajectory. That is, 
it is continuous in time and the differential equation of 
this state will represent its evolution. In Rcfs. [10] and 
[11] we showed that by starting with Eq. (4.4), the time 
derivative of this equation gives diffusive non-Markovian 
SSEs. Thus in this paper we will not reproduce these 
derivations, but instead show that by using our velocity 
operator technique we can rederive the actual trajecto- 
ries for {q k (t)}. This shows that diffusive non-Markovian 
SSEs have a modal interpretation. If fact, because the 
orthodox interpretation only gives an interpretation for 
the solutions of non-Markovian SSE at a particular time 
(time of measurement), we believe that the only non- 
trivial interpretation of non-Markovian SSEs is a modal 
interpretation. 
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Before we consider specific unravelings we would like 
to note that the velocity field, Eq. (2.25), can be written 
in terms of the conditioned system state as 



Vk({qk},t)=Re[(^ {qk} (t)\v k ({q k },t)\i> {qk} (t)}], (4.6) 
where 

MUk}, t)\1> {qk} (t)) = ({q k }\v k (t)\V(t))/VN. (4.7) 

This results in the following differential equation for the 
bath values: 



d t q k {t) = Rc[(^ {qk{t)} (t)\v k ({q k (t)},t)\iP {qkm (t))}, 



where \ip {qk (*)}(*)) = \ip{ qk }(t))\ {qk=qk(t)}- 



(4.8) 



A. Position unraveling 

The first unraveling we consider is the position unrav- 
eling. This results when we chose a preferred projective 
measure of the form 

{^fe} = K{x h } = |{Zfc})({zfc}|bath® isys}, (4.9) 

where is the multimodc cigenstate of the position 

operators x k = (aj, + a k )/\/2. 

To simplify the overall equation we define a noise func- 
tion z(t) as 



(4.10) 



Integrating this gives 

x k (t) = x k (0) + f ds[(L) s g k e^ s + (D) S g k e^ s }/V2 
Jo 

(4.13) 

where x k (0) is the random variable chosen from the dis- 
tribution 



p({.T fc },o)=K{.T fc }|{o fc })i 2 =n 



exp(-x^) 



(4.14) 



That is, we have chosen the initial condition ^(t)) = 
|{0 fc })|^(0)). 

To show that Eq. (4.8) does give the same trajecto- 
ries for the values {x k (t)} as in Eq. (4.13), we apply 
the Hamiltonians defined in Sec. Ill to Eq. (2.26), with 
q k = x k ■ This gives 



v k (t) = [g^L + gke-^tfyV 7 !, 



(4.15) 



as [x k ,d k ] = — 1/\/2 and [xfc,a^] = 1/V%- Substituting 
this into Eq. (4.6) gives a velocity field of the form 



v k ({x k },t) = [g^^^L^^t)) + g k 

x(i> {Xk} (t)\fr\iP {Xk} (t))]/V2. (4.16) 

Thus Eq. (4.8) immediately reproduces Eq. (4.12), 
thereby confirming that the modal theory does give the 
same non-Markovian SSEs, as found with the orthodox 
theory 



B. Quadrature unraveling 



In Rcf. [11] we showed that the non-Markovian SSE for 
the position unraveling is 

d t \Mt)) = {-~H iat (t) + (L-(L) t )z*(t)-(L 

-(L) t )[B z (t) + D z (t)} + {(L - (L) t )[B z {t) 
+t) z {t)}) t -{D -{D) t )[A z {t)+C z {t)] 
+<(Zt -(tf) t )[A z (t) + C z (t)]) t 

x}|V»(*)>, (4-11) 

where |^(*)) = and (L) t = z {t)\L\^ z (t)) . 

The four operators A z (t) , B z (t) , C z (t) and D z (t) are de- 
fined as ansatzen to functional derivatives. It turns out 
that in general these operators are not solvable. The per- 
turbation techniques outlined in Refs. [15] and [16] can 
be applied to this non-Markovian SSE to give a pertur- 
bative solution. Given that there is no Markovian limit 
to this equation, however, it is unclear whether such per- 
turbativc methods would be effective. 

We also showed, after considerable effort, that 

d t x k (t) = [(L^gke* 1 ** + {L^ t g k e- inkt ]/V2. (4.12) 



The next unraveling we consider is what we call the 
quadrature unraveling. In Ref. [10] we show that this 
unraveling only exists for certain environments, such that 
for every mode k in the bath there exists another mode, 
which we can label —k, such that Q- k = —£l k and 
g_fc = g k . These simply imply that the modes coupled 
to the system come in symmetric pairs about the sys- 
tem frequency Q. The form of the preferred projective 
measure for this unraveling is 

= *{X+,Y-} = l{^ + ,^})({^fe + ,^r}lbath®i sys }, 

(4.17) 

where \{X+,Y k ~}) = n fc >o \K ,Y k ~) where \X+,Y k 7) is 
the two mode entangled (EPR) state 



dx' 



-k 



V2 



(4.18) 

Here \{X k + x')/\/2} k is an eigenstate of x k , and simi- 
larly \(X£ -x')/V2)-k for x- k - Equation (4.18) is an 
eigenstate of both the operators 



X 



£ = (x k + X 



(Vk 



k 

y-k 



2. 

2, 



(4.19) 
(4.20) 
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where x k and y k are the quadratures of d k : 
afe = (ifc + iyk)/V2. 



(4.21) 



As in the position unraveling we define a noise function 
z(t) as 

z(t) =^2 5fc [X+(t)cos(Q fc t)+F-(i)sin(Q fe t)], (4.22) 

fc>0 

which by definition is real. In Rcf. [10] we showed that 
the non-Markovian SSE for the quadrature unraveling is 



-j-H int {t) - (L x - (L x ) t )Q z (t) 



+ ((L X - (L x ) t )Q z (t) 



+z(t)(L-(L) t ) \Mt)), 



(4.23) 



L + Lt. 



where \ip z (t)) = ^{x+^.Y' (t)}i 1 )) and L * 
Again the operator Q z (t) is an ansatz to a functional 
derivative. In Ref. [16] we outlined a perturbation tech- 
nique for finding this operator, if an exact solution cannot 
be found. 

The differential equations for X k (t) and Y k ~(t) were 
shown, using the method of characteristics, to be 

d t X+{t) = g k cos{n k t)(L x ) t , (4.24) 
d t Y~{t) = g k sm(n k t){L x ) t . (4.25) 

Integrating these differential equation from time to t 
we get 



X+(t) = X+(0)+ / g k cos(n kS )(L x ) s ds, (4.26) 
Jq 

**"(*) = Y k ~(0)+ f g k wn(Sl k 8){Lt) a d8. (4.27) 



As in the position unraveling the random variables 
A^O) and ^"(O) are chosen from the initial distri- 
bution. For this unraveling and the initial condition 
!*(£)) = |{0 fc })|V>(0)) this distribution is 



p({x+, r fe -},o) = n 



(4.28) 



fe>0 



To show that Eqs. (4.24) and (4.25) can be derived 
from the modal theory (velocity operator technique) we 
apply the Hamiltonians defined in Sec. Ill to Eq. (2.26). 
For this unraveling the set of velocity operators will 
be the union of {v k } and {v k }, where 

«*(*) = ^[^fcAni] =g k L m cos(fl k t), (4.29) 



**(*) 



g k L x sm(Sl k t), (4.30) 



which arc both real by definition. Substituting these ve- 
locity operators into Eq. (4.6) gives 

4({ X k> Y k},t) = 9k{^ {x+ Y - } (t)\L x \4, {x+ Y - } {t)) 
x cos(rj fe t), (4.31) 

4({X^Y k },t) = g k ^ {x + Y - } (t)\L x \^ {x + Y - } (t)) 
x sin(fifct). (4.32) 

Thus Eq. (4.8) simply yields Eqs. (4.24) and (4.25). Thus 
the modal theory gives the correct non-Markovian SSE. 

C. Coherent unraveling 

The last unraveling we consider is the coherent state 
unraveling. This non-Markovian SSE was first presented 
by Diosi, Gisin, and Strunz [5]. In Ref. [10] we showed 
that it could be derived in the orthodox interpretation by 
considering a bath measurement in terms of the Husimi 
POM [35]. This has POM elements 



1 



F M = — \{ a k})({ak}l 



(4.33) 



where a k \a k ) = a k \a k ) and K is the total number of 
modes. The noise function z(t) for this unraveling is 



z(t) = ^2g k a k (t)c 



(4.34) 



In Refs. [5, 10] it is shown that the non-Markovian SSE 
for the coherent unraveling is 



dt\Mt)) = 



[{ti-(L%)C z {t)) t 
+z*{t)(L-(L) t )]\^ z {t)) : 



(4.35) 



where \ip z (t)) = \^{ ak (t)}(t)) ■ The operator C z {t) as 
in the previous unravelings is an ansatz to a functional 
derivative. As in the quadrature unraveling perturbation 
techniques exist for finding a perturbative solution for 
this ansatz [15, 16]. 

Using the same procedure as the other unravelings, the 
differential equation for a k (i) is 



dta* k (t) = g k e 
which integrates to give 



a* k {t) = 4(0) + / g k e 
Jo 



(4.36) 



(4.37) 



For an initial vacuum bath state, the random variable 
a£(0) is defined by the initial distribution 



p({a k }, o)=n 



-Kl 2 



(4.38) 
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To show that the modal theory can be used to describe 
this non-Markovian SSE we have to find the projector in 
K. which is equivalent to the POM elements defined in 
Eq. (4.33). In Ref. [32] we showed that for a single mode 
this projector is \x + 1 y~)(x + ,y~\ where 



/2w 



V2 



x + + x 1 



V2 



^ (4.39) 

where the states in the integrand are x states. Thus the 
multimode projector used to define this unraveling is 

tife} = n {a&} = \{x+,y k })({x+,y k }\ blith+aux ®i sys , 

(4.40) 

where a k is defined by 

a k = xt + iy^. (4.41) 
This allows us to define the operator A k , such that 

A k \x + ,y-) = a k \x+,y-}, (4.42) 



A k = x k +iy k = a k + b\, 



(4.43) 



which is a normal operator [36]. Here x k and y k are 
defined as 



i+ = [a k + a[ + b k + bi}/2, 



*t _i_ t. i hU /<: (4.44) 
y k = [-ia k + ial + ib k -ib k ]/2, (4.45) 



where b k and b\ are annihilation and creation operators 
which act in 7i a ux- In this enlarged Hilbert space the 
velocity operators are 



vt(t) = -tWI J-^uni ® Lux] 



= [g k e in * t L + g k e- in - t D]/2 1 (4.46) 

% * >. 

V k (t) = ~j\Vk i H uni® laux] 

= [-13,6^1 + 19,6-^1^/2. (4.47) 
With these velocity operators the velocity fields become 

<(K.»fc}»*) = ^{^^mgke^L + g.e-^ 
xDM {xty - } (t))/2 7 (4.48) 

«fc({4»»fc}.*) = ^ {x +, y - } (t)\{-me inkt L + tg k 

xe~^D]\i > + (t))/2. (4.49) 

Substituting these into Eq. (4.8) gives 

dtx+(t) = [(L) t9k e inkt + {D) t g k e- inkt ]/2, (4.50) 
dtV k {t) = hmtg^ +i{D) t9k e^ t ]/2. 

(4.51) 

Since a k = x k + iy k we once again easily obtain 
Eq. (4.36) as found from the orthodox theory. 



V. DISCUSSION AND CONCLUSIONS 

Given the success of continuous quantum measurement 
theory [23] in giving a nontrivial interpretation of Marko- 
vian SSEs, it is natural to seek a similar interpretation for 
diffusive non-Markovian SSEs. It turns out that to give a 
meaning to such non-Markovian SSEs we have to consider 
the modal interpretation of quantum mechanics. This is 
because under the orthodox interpretation of quantum 
mechanics, only the solution of the non-Markovian SSE 
at a time t can be given a meaning [10]. It corresponds 
to the state of the system at that time given a measure- 
ment on the bath at that time yielding a particular re- 
sult. The bath cannot be measured continuously because 
non-Markovian systems have a memory, so past measure- 
ments of the bath will in general have a back-action which 
disrupts the system's average evolution. Thus the solu- 
tion at a particular time may have an interpretation, but 
the linking of these solutions at different times does not. 
The trajectory generated by the non-Markovian SSE can 
be regarded only as a numerical tool for calculating a 
conditioned state at a particular time. 

However, under the modal interpretation, in particular 
a view which is closest in line with Bell's beable theory 
[26, 30], we find that non-Markovian SSEs do have a 
nontrivial interpretation. In this interpretation, the bath 
has definite properties even if it is not measured, so the 
backaction problem disappears. The system is treated 
as a purely quantum system which, however, depends 
upon the values of the bath properties: the values of 
the bath hidden variables. The evolution of the system 
is generated by the non-Markovian SSE and the system 
state guides the values of the hidden variables for the 
bath. The bath hidden variables are similar to Bohmian 
hidden variables [33], as they obey a deterministic dif- 
ferential equation with stochastic initial conditions. In 
fact, one unraveling we have considered corresponds to 
Bohmian mechanics of the bath in an interaction frame. 

In this paper we also considered the quadrature and 
coherent unravelings [10]. The different unravelings arc 
determined by choosing which bath observable is to be 
given property status. The noise z(t) which appears in 
the non-Markovian SSE is simply a linear combination of 
the time-dependent values of the bath hidden variables. 
The quadrature unraveling is defined such that z(t) is 
real, while for the coherent unraveling it is complex. In 
the Markovian limit the former becomes the quantum 
trajectory for homodyne detection and the latter be- 
comes the quantum trajectory for heterodyne detection. 
Thus quantum trajectories have both the standard con- 
tinuous measurement interpretation and the above modal 
interpretation. 

In conclusion it seems that to give diffusive non- 
Markovian SSE a satisfying interpretation we must give 
up the orthodox interpretation of quantum mechanics 
and consider the lesser known but equally valid modal 
interpretation. In Refs. [37, 38] it is claimed that the for- 
malism presented there does give a satisfying interpreta- 
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tion of non-Markovian SSEs. If this is correct, we must 
conclude that said formalism is a hidden- variable theory. 
The clarification of this issue is beyond the scope of the 
current paper. Another problem for future work would 
be to determine what choices for objective properties of 
the bath will give rise to SSEs with a Markovian limit. 
Finally, it would be interesting to use the modal theory 
to develop discontinuous non-Markovian SSEs, such as 
those which corresponds to a number state decomposi- 



tion. 
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